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I.  INTRODUCTION 


1.1.  Soft  Active  Materials  for  Soft  Machines 

The  convergence  of  parts  of  biology  and  engineering  has  created  exciting  opportunities 
of  discovery,  invention  and  commercialization.  The  overarching  themes  include  using 
engineering  methods  to  advance  biology,  combining  biology  and  engineering  to  invent  medical 
procedures,  and  mimicking  biology  to  create  engineering  devices. 

Machines  in  engineering  use  mostly  hard  materials,  while  machines  in  nature  are  often 
soft.  What  does  softness  impart  to  the  life  of  animals  and  plants?  A  conspicuous  feature  of  life 
is  to  receive  and  process  information  from  the  environment,  and  then  move.  The  movements 
are  responsible  for  diverse  functions,  far  beyond  the  function  of  going  from  place  to  place.  For 
example,  an  octopus  can  change  its  color  at  an  astonishing  speed,  for  camouflage  and  signaling. 
This  rapid  change  in  color  is  mediated  by  thousands  of  pigment-containing  sacs.  Attach  to  the 
periphery  of  each  sac  are  dozens  of  radial  muscles.  By  contracting  or  relaxing  the  muscles,  the 
sac  increases  or  decreases  in  area  in  less  than  a  second.  An  expanded  sac  may  be  up  to  about  l 
mm  in  diameter,  showing  the  color.  A  retracted  sac  may  be  down  to  about  o.i  mm  in  diameter, 
barely  visible  to  the  naked  eye  M. 

As  another  example,  in  response  to  a  change  in  the  concentration  of  salt,  a  plant  can 
change  the  rate  of  water  flowing  through  the  xylem.  This  regulation  of  flow  is  thought  to  be 
mediated  by  pectins,  polysaccharides  that  are  used  to  make  jellies  and  jams.  Pectins  are  long 
polymers,  crosslinked  into  a  network.  The  network  can  imbibe  a  large  amount  of  water  and 
swell  many  times  its  own  volume,  resulting  in  a  hydrogel.  The  amount  of  swelling  changes  in 
response  to  a  change  in  the  concentration  of  slat.  The  change  in  the  volume  of  the  hydrogel 
alters  the  size  of  the  microchannels  in  the  xylem,  regulating  the  rate  of  flow[2]. 

The  above  examples  are  concerned  with  animals  and  plants.  But  many  more  examples 
are  everywhere  around  and  inside  us.  Consider  the  accommodation  of  the  eye,  the  beating  of 
the  heart,  the  sound  shaped  by  the  vocal  folds,  and  the  sound  in  the  ear.  Abstracting  these 
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biological  soft  machines,  we  may  say  that  a  stimulus  causes  a  material  to  deform,  and  the 
deformation  provides  a  function  (Fig.  l).  Connecting  the  stimulus  and  the  function  is  the 
material  capable  of  large  deformation  in  response  to  a  stimulus.  We  call  such  a  material  a  soft 
active  material  (SAM). 

An  exciting  field  of  engineering  is  emerging  that  uses  soft  active  materials  to  create  soft 
machines.  Soft  active  materials  in  engineering  are  indeed  apt  in  mimicking  the  salient  feature 
of  life:  movements  in  response  to  stimuli.  An  electric  field  can  cause  an  elastomer  to  stretch 
several  times  its  length.  A  change  in  pH  can  cause  a  hydrogel  to  swell  many  times  its  volume. 
These  soft  active  materials  are  being  developed  for  diverse  applications,  including  soft  robots, 
adaptive  optics,  self-regulating  fluidics,  programmable  haptic  surfaces,  and  oilfield 
management13-8!. 

Research  in  soft  active  materials  has  once  again  brought  mechanics  to  the  forefront  of 
human  creativity.  The  familiar  language  finds  new  expressions,  and  deep  thoughts  are 
stimulated  by  new  experience.  To  participate  in  advancing  the  field  of  soft  active  materials  and 
soft  machines  effectively,  mechanicians  must  retool  our  laboratories  and  our  software,  as  well  as 
adapt  our  theories. 

The  biological  phenomena,  as  well  as  the  tantalizing  engineering  applications,  have 
motivated  the  development  of  theories  of  diverse  soft  active  materials,  including  dielectric 
elastomers^'31,  elastomeric  gelsll4_19],  polyelectrolytes[20>21],  pH-sensitive  hydrogels122-24!,  and 
temperature  sensitive  hydrogels [25].  The  theories  attempt  to  answer  commonly  asked  questions. 
How  do  mechanics,  chemistry,  and  electricity  work  together  to  generate  large  deformation? 
What  characteristics  of  the  materials  optimize  their  functions?  How  do  molecular  processes 
affect  macroscopic  behavior?  How  efficiently  can  a  material  convert  energy  from  one  form  to 
another?  The  theories  are  being  implemented  in  software,  so  that  they  can  become  broadly 
useful  in  the  creation  of  materials  and  devices. 


October  25,  2010 


3 


1.2.  Dielectric  Elastomers 


This  review  will  focus  on  one  particular  class  of  soft  active  materials:  dielectric 
elastomers.  All  materials  contain  electrons  and  ions— charged  particles  that  move  in  response 
to  an  applied  voltage.  In  a  conductor,  electrons  or  ions  can  move  over  a  macroscopic  distance. 
By  contrast,  in  a  dielectric,  the  charged  particles  move  relative  to  one  another  by  short  distances. 
The  two  processes— deformation  and  polarization— are  inherently  coupled. 

Fig.2  illustrates  the  principle  of  operation  of  a  dielectric  elastomer  transducer.  A 
membrane  of  a  dielectric  elastomer  is  sandwiched  between  two  compliant  electrodes.  The 
electrodes  have  negligible  electrical  resistance  and  mechanical  stiffness.  A  commonly  used 
material  for  such  electrodes  is  carbon  crease.  The  dielectric  is  subject  to  forces  and  voltage. 
Charge  flows  through  an  external  conducting  wire  from  one  electrode  to  the  other.  The  charges 
of  the  opposite  signs  on  the  two  electrodes  cause  the  membrane  to  deform.  It  was  discovered 
that  an  applied  voltage  may  cause  dielectric  elastomers  to  strain  over  ioo%[3].  Because  of  this 
large  strain,  dielectric  elastomers  are  often  called  artificial  muscles.  The  discovery  has  inspired 
intense  development  of  dielectric  elastomers  as  transducers  for  diverse  applications126-281. 

This  review  focuses  on  the  theory  of  dielectric  elastomers.  Section  II  describes  the 
thermodynamics  of  a  transducer  of  two  independent  variations.  Emphasis  is  placed  on  basic 
ideas:  states  of  the  transducer,  cyclic  operation  of  the  transducer,  region  of  allowable  states, 
equations  of  state,  stability  of  a  state,  and  nonconvex  free-energy  function.  These  ideas  are 
described  in  both  analytical  and  geometrical  terms.  Section  III  develops  the  theory  of 
homogeneous  fields.  After  setting  up  a  thermodynamic  framework  for  electromechanical 
coupling,  we  consider  several  specific  material  models:  a  vacuum  as  an  elastic  dielectric  of 
vanishing  rigidity,  incompressible  materials,  ideal  dielectric  elastomers,  and  electrostrictive 
materials.  Section  IV  applies  nonequilibrium  thermodynamics  to  dissipative  processes,  such  as 
viscoelasticity,  dielectric  relaxation,  and  electrical  conduction.  Section  V  discusses 
electromechanical  instability,  both  as  a  mode  of  failure  and  as  a  means  to  achieve  giant  voltage- 
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induced  deformation.  Section  VI  outlines  the  theory  of  inhomogeneous  fields.  A  variational 
statement  is  formulated  as  the  basis  for  the  finite  element  method.  The  associated  partial 
differential  equations  are  summarized. 


II.  THERMODYANMICS  OF  A  TRANSDUCER 
2.1.  States  of  a  Transducer 

Fig.  3  illustrates  a  transducer,  consisting  of  a  dielectric  separating  two  electrodes.  The 
transducer  is  subject  to  a  force  P,  represented  by  a  weight.  The  two  electrodes  are  connected 
through  a  conducting  wire  to  a  voltage  ®,  represented  by  a  battery.  The  weight  moves  by 
distance  /,  and  the  battery  pumps  charge  Q  from  one  electrode  to  the  other. 

The  transducer  is  capable  of  two  independent  variations.  Consequently,  states  of  the 
transducer  can  be  represented  geometrically  on  a  plane.  The  two  coordinates  of  the  plane  may 
be  chosen  from  variables  such  as  P,  ® ,  /  and  Q.  A  point  in  the  plane  represents  a  state  of  the 
transducer.  For  example,  Fig.  4  shows  the  force  P  and  the  displacement  /  as  the  coordinates  of 
the  plane.  Plotted  on  the  plane  are  the  force-displacement  curves  of  the  transducer.  Each  force- 
displacement  curve  is  measured  under  the  condition  that  the  two  electrodes  are  subject  to  a 
constant  voltage  during  deformation. 

If  the  voltage  between  the  electrodes  can  be  varied,  the  force  and  the  displacement  may 
be  changed  independently.  When  the  transducer  is  subject  to  a  constant  weight,  but  the  voltage 
is  changed  from  o  to  ®,  ,  the  transducer  changes  from  state  A  to  state  B,  lifting  the  weight. 
When  the  displacement  /  is  held  constant,  a  change  in  the  voltage  causes  the  transducer  to 
change  from  state  A  to  state  C,  with  accompanying  change  in  the  force. 

We  could  have  also  plotted  on  the  ( l,P )  plane  curves  of  constant  values  of  charge.  Each 
curve  of  a  constant  charge  is  the  force-displacement  curve  of  the  transducer  measured  under  the 
open-circuit  condition,  when  the  two  electrodes  maintain  a  fixed  amount  of  charge. 
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Following  Gibbs’s  graphical  method  for  a  thermodynamic  system  of  fluids[29],  we  may 
choose  any  two  of  the  four  variables  P,  ® ,  /  and  Q  as  coordinates.  Each  choice  represents  the 
transducer  on  a  different  plane.  All  these  planes  represent  the  same  states  of  the  transducer, 
because  the  transducer  is  capable  only  two  independent  variations.  Nonetheless,  different 
planes  emphasize  different  attributes  of  the  states.  For  example,  the  (P,®)  plane  may  be  used 
to  indicate  loading  conditions,  while  the  ( l,Q )  plane  may  be  used  to  indicate  kinematic 
conditions.  The  (Z,®)  plane  is  often  used  to  report  the  voltage -induced  deformation,  while  the 
(P,Q)  plane  may  be  used  to  report  force-induced  charge. 

2.2.  Cyclic  Operations  of  a  Transducer 

Most  applications  involve  cyclic  changes  of  the  state  of  the  transducer.  A  particular  cycle 
of  states  is  illustrated  in  Fig.  5,  on  the  (Q,®)  plane.  The  voltage  and  the  charge  of  the 
transducer  can  be  changed  independently.  In  changing  from  state  A  to  state  B,  the  transducer  is 
connected  to  a  battery  of  a  low  voltage,  ®,  ;  a  change  in  the  applied  force  reduces  the  spacing 
between  the  two  electrodes,  causing  the  charge  on  the  electrodes  to  increase.  In  changing  from 
state  B  to  state  C,  the  transducer  is  under  an  open-circuit  condition  and  the  electrodes  maintain 
the  constant  charge  Qn ;  a  change  in  the  applied  force  increases  the  spacing  between  the  two 
electrodes,  raising  the  voltage  to  ®H .  In  changing  from  state  C  to  state  D,  the  transducer  is 
connected  to  a  battery  of  the  high  voltage  ®H ;  a  change  in  the  applied  force  increases  the 
spacing  between  the  two  electrodes,  causing  the  charge  on  the  two  electrodes  to  decrease.  In 
changing  from  state  B  to  state  C,  the  transducer  is  under  an  open-circuit  condition  and  the 
electrodes  maintain  the  constant  charge  QL  ;  a  change  in  the  applied  force  decreases  the  spacing 
between  the  two  electrodes,  lowering  the  voltage  to  ® L . 

The  cycle  receives  mechanical  work,  and  pumps  charge  from  a  low  voltage  to  a  high 
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voltage.  Such  a  cycle  describes  a  generator,  harvesting  electric  energy  by  receiving  mechanical 
work  from  the  environment,  such  as  the  work  done  by  an  animal  or  human  during  walking,  and 
the  work  done  by  ocean  waves. 

Indeed,  a  closed  curve  of  any  shape  on  the  (Q,®)  plane  represents  a  cyclic  operation  of 
the  transducer.  The  amount  of  energy  converted  per  cycle  is  given  by  the  area  of  the  cycle  on 
the  (Q,®)  plane.  When  the  states  cycle  counterclockwise  on  the  (Q,®)  plane,  the  transducer  is 
a  generator,  converting  mechanical  energy  to  electrical  energy.  When  the  states  cycle  clockwise 
on  the  (Q,®)  plane,  the  transducer  is  an  actuator,  converting  electrical  energy  to  mechanical 
energy.  The  (Z,P)  plane  can  also  be  used  to  evaluate  energy  of  conversion  per  cycle. 

2.3.  Modes  of  Failure  and  Region  of  Allowable  States 

A  transducer  may  fail  in  multiple  modes,  such  as  mechanical  rupture,  electrical 
breakdown,  electromechanical  instability,  and  loss  of  tension[3<>32f  The  critical  condition  for 
each  mode  of  failure  can  be  represented  on  the  (Q,®)  plane  by  a  curve.  Curves  of  all  modes  of 
failure  bound  in  the  plane  a  region,  which  we  call  the  region  of  allowable  states  of  the  transducer. 
Such  graphic  methods  have  been  used  to  optimize  actuators  [33’34]  and  calculate  the  maximal 
energy  of  conversion  for  generators [35'37].  Fig.  6  shows  an  example[35]. 

2.4.  Equations  of  State 

On  dropping  a  small  distance  SI ,  the  weight  does  work  Pd  .  On  pumping  a  small 
amount  of  charge  SQ  ,  the  battery  does  work  <£>SQ  .  The  force  is  work-conjugate  to  the 
displacement,  and  the  voltage  is  work-conjugate  to  the  charge.  We  will  analyze  isothermal 
processes,  and  remove  temperature  from  explicit  consideration.  Denote  the  Helmholtz  free 
energy  of  the  transducer  by  F.  When  the  transducer  equilibrates  under  the  applied  force  and 
the  applied  voltage,  the  change  in  the  free  energy  of  the  transducer  equals  the  sum  of  the  work 
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done  by  the  weight  and  the  work  done  by  the  battery: 

SF  =  PSl+®SQ.  (l) 

This  condition  of  equilibrium  holds  for  arbitrary  small  variations  SI  and  SQ  .  The  displacement 
and  the  charge  are  independent  variables. 

The  two  independent  variables  ( l,Q )  characterize  the  state  of  the  transducer.  The 
Helmholtz  free  energy  of  the  transducer  is  a  function  of  the  two  independent  variables: 

F  =  F{l,Q).  (2) 


Associated  with  small  variations  SI  and  SQ  ,  the  free  energy  varies  by 

81  8Q 


A  comparison  of  (l)  and  (3)  gives 


\8F(l,Q)  p. 

SL  + 

% 

O 

1 

© 

_ 1 

81  \ 

L  SQ  J 

SQ  =  o. 


(3) 


(4) 


When  the  transducer  equilibrates  with  the  weight  and  the  battery,  the  condition  of  equilibrium 
(4)  holds  for  independent  and  arbitrary  variations  SI  and  SQ .  Consequently,  in  equilibrium, 
the  coefficients  of  the  two  variations  in  (4)  both  vanish,  giving 

pFim,  (5) 

81 

O  FFLQ).  (6) 

dQ 

Once  the  free-energy  function  F(l,Q)  is  known,  (5)  and  (6)  express  P  and  ®  as  functions  of  / 
and  Q.  That  is,  the  two  equations  give  the  force  and  voltage  needed  to  cause  a  certain 
displacement  and  a  certain  charge.  The  two  equations  (5)  and  (6)  constitute  the  equations  of 
state  of  the  transducer. 

Equation  (5)  can  be  used  to  determine  the  free-energy  function  from  the  force- 
displacement  curves  of  the  transducer  measured  under  the  open-circuit  conditions,  when  the 
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electrodes  maintain  constant  charges.  For  each  value  of  Q,  the  free  energy  is  the  area  under  the 
force-displacement  curve.  Similarly,  (6)  can  be  used  to  determine  the  free-energy  function  from 
the  voltage-charge  curves  of  the  transducer.  As  mentioned  before,  (l,P)  and  (Q,®)  are 
convenient  planes  to  represent  the  states  of  the  transducer  when  we  wish  to  highlight  work  and 
energy. 

As  an  illustration,  consider  a  parallel-plate  capacitor— two  plates  of  electrodes  separated 
by  a  thin  layer  of  a  vacuum  (Fig.  7).  The  separation  /  between  the  two  electrodes  may  vary,  but 
the  area  A  of  either  electrode  remains  fixed.  Recall  the  elementary  fact  that  the  amount  of 
charge  on  either  electrode  is  linear  in  the  voltage: 


0  = 


(7) 


where  s0  the  permittivity  of  the  vacuum.  Inserting  (7)  into  (6)  and  integrating,  we  obtain  the 
free-energy  function 


F(l.Qh 


IQ2 

2£qA' 


(8) 


Inserting  (8)  into  (5),  we  obtain  that 


(9) 


Equations  (7)  and  (9)  constitute  the  equations  of  state  of  the  parallel-plate  capacitor.  They  are 
readily  interpreted.  The  applied  voltage  causes  charge  to  flow  from  one  electrode  to  the  other, 
so  that  one  electrode  is  positively  changed,  and  the  other  is  negatively  charged.  Equation  (7) 
relates  the  charge  to  the  applied  voltage.  The  oppositely  charged  electrodes  attract  each  other. 
To  maintain  equilibrium,  a  force  need  be  applied  to  the  electrode.  Equation  (9)  relates  the 
applied  force  to  the  charge. 

Define  the  electric  field  by  E  =  ®/Z  and  the  stress  by  a-P/A.  Rewrite  (9)  as 
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(10) 


E2. 


This  equation  gives  the  stress  needed  to  be  applied  to  the  electrodes  to  counteract  the 
electrostatic  attraction.  This  stress  is  known  as  the  Maxwell  stress. 


2.5.  Stability  of  a  State 

The  equations  of  state,  (5)  and  (6),  are  in  general  nonlinear.  If  the  transducer  operates 
in  the  neighborhood  of  a  particular  state  ( l,Q ),  the  equations  of  state  can  be  linearized  and 
written  in  an  incremental  form: 


SP  = 


S®  = 


dl2  dldQ 

(n) 

Mf.Q)  ri%Q) 
dQdl  dQ2 

(12) 

We  call  d2F(l,Q)/  dl2  the  mechanical  tangent  stiffness  of  the  transducer,  and  d2F(l,Q)/  dQ2  the 
electrical  tangent  stiffness  of  the  transducer.  The  two  electromechanical  coupling  effects  are 
both  characterized  by  the  same  cross  derivative,  d2F(l,Q)/  dldQ  =  d2F(l,Q)/  dQdl .  Consequently, 


the  two  electromechanical  coupling  effects  reciprocate.  The  matrix 


H(Z,Q) 


d*F{l,Q)  d2F(l,Q) 
8l2  dldQ 

d2F{l,Q )  d2F{l,Q) 

dldQ  dQ2 


(13) 


is  known  as  the  Hessian  of  the  free-energy  function  F(l,Q). 

As  mentioned  above,  a  state  of  the  transducer  can  be  represented  by  a  point  in  the  ( l,Q ), 
as  well  as  by  a  point  in  the  (P,®)  plane.  For  the  same  state  of  the  transducer,  the  point  in  the 
(Z,Q)  plane  is  mapped  to  the  point  in  the  (P,®)  plane  by  the  equations  of  state,  (5)  and  (6).  The 
mapping  may  not  always  be  invertible.  That  is,  given  a  pair  of  the  loads  (P,®),  the  equations  of 
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state  may  not  be  invertible  to  determine  a  state  ( l,Q ).  For  example,  (n)  and  (12)  are  not 
invertible  when  the  Hessian  is  a  singular  matrix,  detH  =  o .  This  singularity  may  be  understood 
in  terms  of  thermodynamics. 

The  transducer  and  the  loading  mechanisms  (i.e.,  the  weight  and  the  battery)  together 
constitute  a  thermodynamic  system.  The  free  energy  of  the  system  is  the  sum  of  the  free 
energies  of  the  individual  parts— the  transducer,  the  weight,  and  the  battery.  The  free  energy 
(i.e.,  the  potential  energy)  of  a  constant  weight  is  -  PI .  The  free  energy  of  a  battery  of  a  constant 
voltage  is  -  OQ  .  Consequently,  the  free  energy  of  the  thermodynamic  system  combining  the 
transducer  and  the  loading  mechanisms  is 

G(l,Q)  =  F(l,Q)-Pl-(t>Q.  (14) 

The  system  has  two  independent  variables,  /  and  Q. 

Thermodynamics  requires  that  the  system  should  reach  a  stable  state  of  equilibrium 
when  the  free-energy  function  G(Z,Q)  is  a  minimum  against  small  changes  in  /  and  Q.  When  the 
weight  moves  by  SI  and  the  battery  pumps  charges  SQ  ,  the  free  energy  of  the  system  varies  by 


5G  = 


wo?)  P 

Sl  + 

O 

1 

© 

_ 1 

dl 

L  8Q  J 

SQ 


(15) 


25/ 2  v  '  dldQ  v"/v'^/  20Q2 

We  have  expanded  the  Taylor  series  up  to  terms  quadratic  in  SI  and  SQ  .  In  a  state  of 
equilibrium,  the  coefficients  of  the  first-order  variations  vanish,  recovering  the  equations  of 
state  (5)  and  (6).  To  ensure  that  this  state  of  equilibrium  minimizes  G,  the  sum  of  the  second- 
order  variations  must  be  positive  for  arbitrary  combination  of  SI  and  SQ .  That  is,  a  state  of 
equilibrium  is  stable  if  the  Hessian  of  the  free  energy  of  the  transducer,  H (l,Q),  is  positive- 
definite.  The  two-by-two  matrix  is  positive-definite  if  and  only  if 


wtw)>  w(w)>0, 


5/2 


5Q2 


" d2F{l,Q)~ 

> 

~  dzF(l,Q) 

5/2 

5Q2 

dldQ 

(16) 
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When  the  Hessian  of  the  free  energy  function  is  positive-definite,  the  function  F(l,Q)  is  convex 
at  this  state  ( l,Q ). 

As  an  illustration,  consider  the  parallel-plate  capacitor  again.  Given  the  free-energy 
function  (8),  the  second  derivatives  are 

d2F{l,Q)  d*F{l,Q)  l  \d*F{l,Q)]  Q 
dF  ’  8Q2  s0A  ’  [  818Q  \  saA ' 

Consequently,  the  Hessian  is  not  positive-definite  at  any  state  of  equilibrium.  That  is,  the 

parallel-plate  capacitor  subject  to  a  constant  force  and  a  constant  voltage  cannot  reach  a  stable 

state  of  equilibrium.  The  conclusion  is  readily  understood.  The  weight  is  independent  of  the 

separation  between  the  plates,  but  the  electrostatic  attractive  force  increases  as  the  separation 

decreases.  Subject  to  a  fixed  weight,  the  two  plates  will  be  pulled  apart  if  the  voltage  is  small, 

and  will  be  pulled  together  if  the  voltage  is  small. 

The  capacitor  can  be  stabilized  by  a  modification  of  the  loading  mechanisms.  For 

example,  we  can  replace  the  weight  with  a  spring  that  restrains  the  relative  movement  of  the 

plates.  Let  Kbe  the  stiffness  of  the  spring,  and  lQ  be  the  separation  between  the  electrodes  when 

the  spring  is  unstretched,  so  that  the  force  in  the  spring  is  P  -  K(l0  -Z).  The  free  energy  of  the 

system  is  the  sum  of  the  free  energies  of  the  capacitor,  the  spring  and  the  battery: 

G(l,Q) - —  K{1  - 1  )2  -OQ  .  (18) 

2  SaA  2 

In  a  state  of  equilibrium,  the  first  derivatives  of  G(l,Q)  vanish,  giving  the  same  equations  of 
state  as  (to)  and  (12).  The  state  of  equilibrium  is  stable  if  and  only  if  the  Hessian  of  G(l,Q)  is 
positive-definite.  The  second  derivatives  of  the  function  G(l,Q )  are 

d*G{l,Q)_K  82G{l,Q)_  l  d2G{l,Q)_  Q 

8l2  ’  8Q2  s0A  ’  818Q  s0A  '  L 

A  state  of  equilibrium  (l,Q)  is  stable  if  and  only  if 
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Kl 

8,A 


f  Q_" 
\e*aj 


(20) 


Thus,  the  transducer  is  stable  when  the  applied  voltage  is  sufficiently  small. 


2.6.  Nonconvex  Free-Energy  Surface 

Following  Gibbs[38],  we  may  interpret  above  analytical  statements  geometrically. 
Consider  a  three-dimensional  space  with  ( l,Q )  as  the  horizontal  plane,  and  F  as  the  vertical  axis. 
In  this  space,  the  Helmholtz  free  energy  F(l,Q)  is  represented  by  a  surface.  A  pair  of  the  loads 
(P,  o)  is  represented  by  an  inclined  plane  passing  through  the  origin  of  the  space,  with  P  as  the 
slope  with  respect  to  the  /  axis,  and  ®  the  slope  of  the  tangent  plane  with  respect  to  the  Q  axis. 
By  definition  (14),  the  function  G(l,Q)  is  the  vertical  distance  between  the  Helmholtz  free 
energy  surface  and  the  inclined  plane.  Thermodynamics  dictates  that  this  vertical  distance 
G(l,Q )  should  minimize  when  the  transducer  equilibrates  with  the  loads. 

When  the  loads  (P,®)  are  given,  we  may  picture  a  plane  simultaneously  parallel  to  the 
inclined  plane  and  tangent  to  the  free-energy  surface  F(l,Q).  That  is,  the  slope  of  the  tangent 
plane  with  respect  to  the  /  axis  is  P,  and  the  slope  of  the  tangent  plane  with  respect  to  the  Q  axis 
represents  is  ®.  From  the  geometry,  it  is  evident  that  the  state  (l,Q)  of  the  tangent  point 
minimizes  the  vertical  distance  G(l,Q).  Also  evident  from  the  geometry,  G(l,Q )  is  minimum 
only  if  the  free-energy  surface  F(l,Q )  in  the  neighborhood  of  the  state  (/,£))  is  above  the  tangent 
plane— that  is,  the  surface  F(l,Q )  is  convex  at  the  state  (/,Q). 

When  the  loads  (P,®)  change  gradually,  so  are  the  slopes  of  the  inclined  plane. 
Consequently,  as  the  loads  (P,®)  change  gradually,  the  associated  tangent  plane  rolls  along  the 
free-energy  surface.  If  the  free-energy  surface  F(l,Q)  is  globally  convex,  every  tangent  plane 
touches  the  surface  at  only  one  point,  and  only  one  state  ( l,Q )  is  associated  with  a  pair  of  given 
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loads  (P,®).  By  contrast,  if  part  of  the  free-energy  surface  F(l,Q)  is  concave,  a  tangent  plane 
may  touch  the  surface  at  two  points,  and  two  states  (Z,Q)  are  associated  with  a  pair  of  given 
loads  (P,®). 

It  was  discovered  that  the  free  energy  functions  for  dielectric  elastomers  are  typically 
nonconvex[39].  Associated  with  a  given  set  of  loads,  more  than  one  states  of  equilibrium  may 
exist.  The  practical  significance  of  nonconvex  free-energy  functions  will  be  discussed  later  in 
connection  with  electromechanical  instability. 

III.  HOMOGENEOUS  FIELD 

We  now  develop  a  field  theory  of  deformable  dielectrics.  The  field  theory  assumes  that  a 
material  is  a  sum  of  many  small  pieces,  and  the  field  in  each  small  piece  is  homogeneous.  This 
assumption  enables  us  to  define  quantities  per  unit  length,  per  unit  area,  and  per  unit  volume. 
This  section  focuses  on  the  homogeneous  field  of  a  small  piece,  and  Section  VI  considers 
inhomogeneous  field  of  a  body  by  summing  up  small  pieces. 

This  section  begins  by  setting  up  a  thermodynamic  framework  for  electromechanical 
coupling.  We  then  consider  several  specific  material  models:  a  vacuum  as  an  elastic  dielectric  of 
vanishing  rigidity,  incompressible  materials,  ideal  dielectric  elastomers,  and  electrostrictive 
materials. 

3.1.  Equations  of  State 

With  reference  to  Fig.  2,  consider  a  block  of  an  elastic  dielectric,  sandwiched  between 
two  compliant  electrodes.  In  the  reference  state,  the  dielectric  is  subject  to  no  forces  and 
voltage,  and  is  of  dimensions  L, ,  L2  and  P3.  In  the  current  state,  the  dielectric  is  subject  to 

forces  P1 ,  P,  and  P3 ,  and  the  two  electrodes  are  connected  to  a  battery  of  voltage  ®  through  a 
conducting  wire.  In  the  current  state,  the  dimensions  of  the  dielectric  become  /, ,  l2  and  l3 ,  the 
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two  electrodes  accumulate  electric  charges  ±Q  ,  and  the  Helmholtz  free  energy  of  the 
membrane  is  F . 

When  the  dimensions  of  the  dielectric  change  by  d/, ,  Sl2  and  Sl3,  the  forces  do  work 
PA+P&  +  P3&3 .  When  a  small  quantity  of  charge  SQ  flows  through  the  conducting  wire,  the 
voltage  does  work  OdQ  •  When  the  dielectric  equilibrates  with  the  forces  and  the  voltage,  the 
increase  in  the  free  energy  equals  the  work  done: 

SF  =  PlSll  +  P,Ac%  +  P3Sl3  +  OPQ .  (21) 

The  condition  of  equilibrium  (21)  holds  for  arbitrary  small  variations  of  the  four  independent 
variables,  lt ,  l2 ,  /,,  and  Q. 

Define  stretches  by  k1=l1/L1  ,  T2=Z2/L2  and  T3=Z3/L3  ,  nominal  stresses  by 
s1=P1/(L2L3 )  ,  s2=P2/(L1L3)  and  s3  =  P3/(L1L2)  ,  nominal  electric  field  by  P  =  ®/L3, 

nominal  electric  displacement  by  D  =  Q/(L1L2) ,  and  nominal  density  of  the  Helmholtz  free 
energy  by  W  =  f7/(LlL2L3) .  Also  define  true  stresses  by  cr  =  P  /(ij.)  ,  <j2  =P2/(Z1Z3)  and 
cr3  =P3 /(/,/;,),  true  electric  field  by  E  =  O / /,, ,  and  true  electric  displacement  by  D  =  Q/(lllj. 

The  condition  of  equilibrium  (21)  holds  in  any  current  state.  However,  it  is  convenient 
to  divide  both  sides  of  (21)  by  L^L.^ ,  the  volume  of  the  block  in  the  reference  state.  We  obtain 
that 

SW  —  S1SA1  +  S2SA2  +  S3 SA3  +  ESP)  .  (22) 

The  condition  of  equilibrium  holds  for  arbitrary  small  variations  of  the  four  independent 
variables,  k, ,  k2,  k3  and  D . 

As  a  material  model,  the  nominal  density  of  the  Helmholtz  free  energy  is  prescribed  as  a 
function  of  the  four  independent  variables: 

w=w(t,,t2,4,d).  (23) 
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Inserting  (23)  into  (22),  we  obtain  that 


(  8W  ' 

- s, 

S\+\ 

[  dw  ^ 
- 

S\  + 

sw  1 
- 

S\  + 1 

fdW 

10 

Is**  ) 

3  I 

1 8D  J 

SD  =  o . 


(24) 


This  condition  of  equilibrium  holds  for  any  small  variations  of  the  four  independent  variables. 
Consequently,  when  the  dielectric  is  in  equilibrium  with  the  applied  forces  and  the  applied 
voltage,  the  coefficient  in  front  of  the  variation  of  each  independent  variable  vanishes,  giving 


S3  = 


E  = 


8\ 

8\ 

M 3 

8w(a±,A2,A 3,d) 

db 


(25) 


(26) 


(27) 


(28) 


Equations  (25)-(28)  constitute  yet  another  representation  of  the  condition  of  equilibrium— they 
are  called  the  equations  of  state.  Once  the  free-energy  function  w(At,A2,A3,bj  is  prescribed  as  a 
material  model,  the  equations  of  state  (25)-(28)  give  the  values  of  the  forces  and  voltage  needed 
to  equilibrate  with  the  dielectric  in  the  state  {a1,a2,a3,d). 

In  the  absence  of  the  applied  forces,  the  stresses  in  the  dielectric  vanish.  The  stresses 
are  zero  even  when  the  voltage  causes  the  dielectric  to  deform.  Thus,  when  the  battery  applies  a 
voltage  to  the  dielectric,  the  positive  charge  on  one  electrode  and  the  negative  charge  on  the 
other  electrode  cause  the  dielectric  to  thin  down.  We  simply  report  what  we  have  observed  in 
this  experiment:  the  voltage  causes  the  dielectric  to  deform.  We  do  not  jump  to  the  conclusion 
that  the  voltage  causes  a  compressive  stress.  In  this  regard,  we  view  the  deformation  caused  by 
the  voltage  in  the  same  way  as  we  view  the  deformation  caused  by  a  change  in  temperature: 
both  are  stress-free  deformation,  so  long  as  the  material  is  unconstrained[12]. 
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The  work  done  by  the  battery,  OdQ ,  can  be  written  as 

0><5Q  =  (l,eHL>L,d)=  LtL2L?Edf) .  (29) 

That  is,  when  a  dielectric  deforms,  the  nominal  electric  field  and  the  nominal  electric 
displacement  are  work-conjugate.  By  contrast,  in  terms  of  the  true  electric  field  and  the  true 
electric  displacement,  the  work  done  by  the  battery  is 

O  SQ  =  (lEjsd.lD)  =  l?EDS(lJ.2)+ltlJ.?E5D .  (30) 

For  a  deformable  dielectric,  o  ,  so  that  the  true  electric  displacement  is  not  work- 

conjugate  to  the  true  electric  field[12]. 


3.2.  Vacuum 

As  an  application  of  the  equations  of  state  (25M28),  consider  a  block  of  a  vacuum.  We 
think  of  the  vacuum  as  an  elastic  dielectric  with  vanishing  rigidity,  undergoing  a  homogenous 
deformation  k, ,  A2,  and  \ .  Recall  an  elementary  fact  that  saEz  /  2  is  the  electrostatic  energy  of 

block  divided  by  the  block  in  the  current  state.  Consequently,  the  nominal  density  of  free  energy 
is 


W  =  -e0E%A2I  . 

2 


(31) 


Recall  that  the  true  electric  displacement  relates  to  the  true  electric  field  as  D  =  s0E  ,  and 


relates  to  the  nominal  electric  displacement  as  D  =  D/(\A2).  We  rewrite  (31)  in  terms  of  the 
stretches  and  the  nominal  electric  displacement: 


W I 


(Ai,A2,\,d)~- 


T>21 


2^oV2 


(32) 


Inserting  (32)  into  (25M28),  we  obtain  that 


s.  =  — 


du3 
2*<AX 


(33) 
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s. 


(34) 


D'1! 


16  XX 


S 


3 


D2 

’ 


(35) 


E  = 


DA„ 


s  XX 


(36) 


Equations  (33)-(36)  can  be  expressed  in  terms  of  the  true  quantities  as 


<t,=—e0E\ 

2 


a2=  —  60E\ 

2 


y-»2 

az  =  ~£oE  > 

D  =  s„E- 


(37) 

(38) 

(39) 

(40) 


Equations  (37M39)  recover  the  stresses  obtained  by  Maxwell [40].  They  are  valid  when  the 
electric  field  is  in  direction  3. 

The  Maxwell  stress  is  a  tensor.  We  have  already  interpreted  the  component  of  the 
Maxwell  stress  in  the  direction  of  the  electric  field.  We  now  look  at  the  two  components  of  the 
Maxwell  stress  transverse  to  the  direction  of  the  electric  field.  Fig. 8  illustrates  a  classic 
experiment  of  a  capacitor,  which  is  partly  in  the  air  and  partly  in  a  dielectric  liquid.  The  applied 
voltage  causes  the  liquid  to  rise  to  a  height  h.  The  height  results  from  the  balance  of  the 
Maxwell  stress  and  the  weight  of  the  liquid.  The  Maxwell  stress  parallel  to  the  electrodes  in  the 
air  is  cra  =  -saE2  /  2  >  where  sa  is  the  permittivity  of  the  air.  The  Maxwell  stress  parallel  to  the 


electrodes  in  the  liquid  is  cr,  =  -s,E2  /  2 ,  where  6,  is  the  permittivity  of  the  liquid.  The  electric 
field  near  the  air/liquid  interface  is  distorted,  so  that  the  above  two  formulas  are  correct  only  at 
some  distance  away  from  the  interface.  Because  st  >  sa ,  the  difference  in  the  Maxwell  stresses 
in  the  two  media  will  draw  the  liquid  up  against  gravity.  Examining  the  free-body  diagram,  and 
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balancing  the  electrostatic  forces  with  the  weight  of  the  liquid,  we  obtain  that 
pgh  =  {sl  -  sa  )E2  /  2 ,  where  pg  is  the  weight  per  unit  volume  of  the  liquid. 


3.3.  Incompressibility 

When  an  elastomer  undergoes  large  deformation,  the  change  in  the  shape  is  typically 
much  larger  than  the  change  in  the  volume.  Consequently,  the  elastomer  is  often  taken  to  be 
incompressible— that  is,  the  volume  of  the  material  remains  unchanged  during  deformation, 
l,l,J3  =  L,L2L:j ,  so  that 

W3  = 1  •  (41) 

This  assumption  of  incompressibility  places  a  constraint  among  the  three  stretches.  We  regard 
\  and  A2  as  independent  variables,  rewrite  (41)  as  A3  =  i/(A1A2),  and  express  8A3  in  terms  of 

SA,  and  SA2 : 


SA 


3 


SA,  SA2 

KK  KK 


(42) 


In  terms  of  the  variations  of  the  independent  variables,  the  condition  of  equilibrium  (22) 
becomes 


dW  = 


Q - — 

v  1  KKj 


SA, 


S - — 

v2  KKj 


SA+ESD. 


(43) 


For  an  incompressible  material,  this  condition  of  equilibrium  holds  for  any  small  variations  of 
the  three  independent  variables,  A, ,  A2  and  D  . 

For  an  incompressible  elastic  dielectric,  the  density  of  the  free  energy  is  a  function  of  the 
three  independent  variables: 


(44) 


Inserting  (44)  into  (43),  we  obtain  that 
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Because  <54, ,  SA2  and  SD  are  independent  variations,  the  condition  of  equilibrium  (45)  is 


equivalent  to  three  equations: 


(46) 

(47) 

(48) 


The  four  equations,  (41)  and  (46M48),  constitute  the  equations  of  state  for  an  incompressible 
dielectric  elastomer. 


3.4.  Ideal  Dielectric  Elastomers 

An  elastomer  is  a  three-dimensional  network  of  long  and  flexible  polymers,  held 
together  by  crosslinks  (Fig.  9).  Each  polymer  chain  consists  of  a  large  number  of  monomers. 
Consequently,  the  crosslinks  have  negligible  effect  on  the  polarization  of  the  monomers:  the 
elastomer  can  polarize  nearly  as  freely  as  a  polymer  melt.  The  permittivity  changes  by  only  a 
few  percent  when  the  area  of  a  membrane  of  an  elastomer  is  stretched  25  times  [4l].  As  an 
idealization,  we  may  assume  that  the  dielectric  behavior  of  an  elastomer  is  exactly  the  same  as 
that  of  a  liquid  polymer,  so  that  the  density  of  free-energy  function  takes  the  form 

W  =  Ws{\,K)+^sE\  (49) 

where  E  is  the  true  electric  field,  s  the  permittivity,  and  Ws  (l, ,  A, )  the  free  energy  associated 
with  the  stretching  of  the  elastomer.  The  material  is  also  taken  to  be  incompressible,  U4  - 1  • 
We  call  this  material  model  the  model  of  ideal  dielectric  elastomers. 
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The  true  electric  displacement  relates  to  the  true  electric  field  by  D-sE,  and  relates  to 
the  nominal  electric  field  as  D  =  D/{\\)  .  Consequently,  in  terms  the  nominal  electric 
displacement,  (49)  becomes 

w{\,\,d)=Ws{\,K)+±-  -fr-  .  (50) 

Inspecting  (49)  and  (50),  we  note  that  the  electromechanical  coupling  in  an  ideal  dielectric 
elastomer  is  purely  a  geometric  effect. 

Inserting  (50)  into  (46M48)  and  taking  the  partial  differentiations,  we  express  the 
results  in  terms  of  the  true  quantities: 


_  _  ,8W,fa, 4) 

ai  a3  ~  A  £h  > 

oA1 

(51) 

,  8WS{' \,\)  772 

^  -’-SE  . 

(52) 

D  =  eE. 

(53) 

The  above  expression  shows  that  a  through-thickness  electric  field  adds  a  compressive  stress  of 
magnitude  sE2  in  the  two  in-plane  directions.  This  set  of  equations  of  state  has  been  used 
almost  exclusively  in  all  analyses  of  dielectric  elastomers.  The  equations  are  usually  justified  in 
terms  of  the  Maxwell  stress[3].  Now  we  have  interpreted  these  equations  using  the  model  of 
ideal  dielectric  elastomers.  That  is,  the  Maxwell  stress  is  valid  when  the  dielectric  behavior  of 
the  material  is  liquid-like,  unaffected  by  deformation. 

Observe  that,  in  (51)  and  (52),  the  magnitude  of  the  voltage-induced  stress  is  twice  the 
magnitude  of  the  Maxwell  stress.  This  apparent  difference  is  readily  understood  (Fig.  10). 
Because  the  elastomer  is  taken  to  be  incompressible,  superposition  of  a  state  of  hydrostatic 
stress  does  not  affect  the  state  of  deformation.  Start  from  the  state  of  triaxial  stresses 
(-sE2 /2,-sE2 /2,  +  aE2 /z),  as  derived  by  Maxwell,  a  superposition  of  a  state  of  hydrostatic 
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stress  {+ sE2  /  2,  + sE2  /  2,  + sE2  / 2)  gives  a  state  of  uniaxial  stress  (o,o,  +  £E2)  ,  and  a 
superposition  of  a  state  of  hydrostatic  stress  (-  sE 2  /  2 ,-sE2  /  2 ,-sE2  /  2)  gives  a  state  of  biaxial 
stress  {-  sE2 sE2 ,  o).  For  an  incompressible  material,  the  three  states  of  stress  illustrated  in 
Fig.  10  cause  the  same  state  of  deformation. 

The  free  energy  due  the  stretching  of  the  elastomer,  Ws  (A[  ,A2),  may  be  selected  from  a 
large  menu  of  well-tested  functions  in  the  theory  of  rubber  elasticity.  For  example,  the  neo- 
Hookean  model  takes  the  form 


Ws=^(a:+4+4-3),  (54) 

where  //  is  the  small-strain  shear  modulus. 

In  an  elastomer,  each  individual  polymer  chain  has  a  finite  contour  length.  When  the 
elastomer  is  subject  no  loads,  the  polymer  chains  are  coiled,  allowing  a  large  number  of 
conformations.  Subjected  to  loads,  the  polymer  chains  become  less  coiled.  As  the  loads 
increase,  the  end-to-end  distance  of  each  polymer  chain  approaches  the  finite  contour  length, 
and  the  elastomer  approaches  a  limiting  stretch.  On  approaching  the  limiting  stretch,  the 
elastomer  stiffens  steeply.  This  effect  is  absent  in  the  neo-Hookean  model,  but  is  represented  by 
the  Arruda  and  Boyce  models]  and  the  Gent  model[44].  The  latter  takes  the  form 


where  //  is  the  small-stress  shear  modulus,  and  Jlim  is  a  constant  related  to  the  limiting  stretch. 
The  stretches  are  restricted  as  o  <  (a2  +  T2  +  A2  -  3)/  Jhm  <  1 .  When  (a2  +  T2  +  T2  -  3)/  Jlim  — >  o ,  the 
Taylor  expansion  of  (55)  is  (54).  That  is,  the  Gent  model  recovers  the  neo-Hookean  model  when 
deformation  is  small  compared  to  the  limiting  stretch.  When  (/l2  +  T2  +  T2  -  3)/  Jlim  — >  1 ,  the  free 
energy  diverges,  and  the  elastomer  approaches  the  limiting  stretch. 
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3-5-  Electrostriction 

The  voltage  may  cause  some  dielectrics  to  become  thinner,  but  other  dielectrics  to 
become  thicker  (Fig.  11).  For  dielectrics  that  are  nonpolar  in  the  absence  of  electric  field,  the 
voltage-induced  deformation  has  been  analyzed  by  invoking  stresses  of  two  origins: 
electrostriction  and  the  Maxwell  stress.  The  electrostriction  results  from  the  effect  of 
deformation  on  permittivity. 

As  a  simplest  model  of  electrostriction,  we  expand  the  free-energy  function  W  into  the 
Taylor  series  in  powers  of  E  up  to  the  quadratic  term.  The  expansion  is  written  in  the  form  of 
(49),  but  now  the  permittivity  is  a  function  of  the  stretches: 

s  =  s(A1,A2).  (56) 


The  same  procedure  as  the  above  now  gives  the  following  equations  of  state[45]: 


oA,  2  oA, 


^  dWs{A  A2) 
oAn 


a2ez  gg(vy 

2  dA2 


(58) 


D  =  s(A1,A2)E  . 


(59) 


The  variation  of  the  permittivity  with  stretches  has  been  observed  experimentally[46]. 
Further  measurements  are  needed  to  ascertain  the  practical  significance  of  electrostriction  in 
dielectric  elastomers. 


IV.  NONEQUILIBRIUM  THERMODYNAMIOCS  OF  DYELECTRIC  ELASTOMERS 

An  elastomer  responds  to  forces  and  voltage  by  time-dependent,  dissipative  processes [47_ 
49].  Viscoelastic  relaxation  may  result  from  slippage  between  long  polymers  and  rotation  of 
joints  between  monomers.  Dielectric  relaxation  may  result  from  distortion  of  electron  clouds 
and  rotation  of  polar  groups.  Conductive  relaxation  may  result  from  migration  of  electrons  and 
ions  through  the  elastomer.  This  section  describes  an  approach  to  construct  models  of 
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dissipative  dielectric  elastomers,  guided  by  nonequilibrium  thermodynamics [50]. 

Thermodynamics  requires  that  the  increase  in  the  free  energy  should  not  exceed  the 
total  work  done,  namely, 

SF  <  PlSll  +  P2d2  +  P3SL  3  +  OdQ .  (60) 

For  the  inequality  to  be  meaningful,  the  small  changes  are  time-directed:  8j  means  the  change 
of  the  quantity/from  one  time  to  a  slightly  later  time. 

Divide  both  sides  of  (60)  by  the  volume  of  the  membrane,  fyfyL,,  ,  and  the 
thermodynamic  inequality  becomes 

8W  <  s1SA1  +  s2SZ2  +  s3SZ3  +  E8D  .  (61) 

As  a  model  of  the  dielectric  elastomer,  the  free-energy  density  is  prescribed  as  a  function: 

W  =  w(\,\,\,D,^2,..).  (62) 


We  characterize  the  state  of  a  dielectric  by  Z2,  X3  and  D  ,  along  with  additional  parameters 
Inspecting  (6i),  we  note  that  At ,  A2,  />,  and  D  are  the  kinematic  parameters  through 
which  the  external  loads  do  work.  By  contrast,  the  additional  parameters  (£, £,,••)  are  not 
associated  with  the  external  loads  in  this  way.  These  additional  parameters  describe  the  degrees 
of  freedom  associated  with  dissipative  processes,  and  are  known  as  internal  variables. 

Inserting  (62)  into  (61),  we  rewrite  the  thermodynamic  inequality  as 


^<0. 


(63) 


As  time  goes  forward,  this  thermodynamic  inequality  holds  for  any  change  in  the  independent 
variables  {z l9>i2,/l3,D,4i,£2,..).  We  next  specify  a  model  consistent  with  this  inequality. 


We  assume  that  the  system  is  in  mechanical  and  electrostatic  equilibrium,  so  that  in  (63) 


the  factors  in  front  of  SZ1 ,  SZ2  and  SD  vanish: 
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„  _8w{ai,A2,A3,D,^1,^2,.., 

-  -  > 

(64) 

„  _dw(\,A2,A3,D,^2,..) 

$2  -  .  ? 

(65) 

dX2 

„  _dw(A1,A2,A3,D,^2,..) 

^0  ^  -  ? 

(66) 

3  sa3 

~_8W(A1,A2,A3,D,^,^2,..) 

(67) 

8D 


Equations  (64M67)  constitute  the  thermodynamic  equations  of  state  of  the  dielectric  elastomer. 

Once  the  elastomer  is  assumed  to  be  in  mechanical  and  electrostatic  equilibrium,  the 
inequality  (63)  becomes 


Z 


8w(ai,A2,A3,D,^1,^2,... 


^,<0. 


(68) 


This  thermodynamic  inequality  may  be  satisfied  by  prescribing  a  suitable  relation  between 
(S^lfS^2,..)  and  (dW  /  8^, 8W  /  d%2,..).  For  example,  one  may  adopt  a  kinetic  model  of  the  type 


d^  =  __yM  dw(A1,A2,A3,D,^2,.) 
dt  ^  tj  8$ 


(69) 


Here  Mtj  is  a  positive-definite  matrix,  which  may  depend  on  the  independent  variables 

(a1,a2,a3,d,^2,  ..). 

To  represent  a  dissipative  dielectric  elastomer  using  the  above  approach,  we  need  to 
specify  a  set  of  internal  variables  (£,£>••)  ,  and  then  specify  the  functions 

w(ai  ,  A2 ,  A3 ,  D,  c, ,  c2  >  •  • )  and  Mij  (fy  ,  A2 ,  A3 ,  D,  c , ,  c, , . . )  •  There  is  considerable  flexibility  in 

choosing  kinetic  models  to  fulfill  the  thermodynamic  inequality  (68). 

Viscoelastic  relaxation  is  commonly  pictured  with  an  array  of  springs  and  dashpots, 
known  as  the  rheological  models;  see  a  recent  example[5l].  Similarly,  dielectric  relaxation  is 
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commonly  pictured  with  models  consisting  of  resistors  and  capacitors.  By  contrast,  electrical 
conduction  involves  the  transport  of  charged  species  over  a  long  distance.  Coupled  large 
deformation  and  transport  of  charged  species  are  significant  in  polyelectrolytes[21],  and  will  not 
be  discussed  here. 


V.  ELECTROMECHANICAL  INSTABILITY 

While  all  dielectrics  deform  under  voltage,  the  amount  of  deformation  differs  markedly 
among  different  materials.  Under  voltage,  piezoelectric  ceramics  attain  strains  of  typically  less 
than  1%.  Glassy  and  semi-crystalline  polymers  can  attain  less  than  10%  [52].  Strains  about  30% 
were  observed  in  some  elastomers [53].  In  the  last  decade,  strains  over  100%  have  been  achieved 
in  several  ways,  by  pre-stretching  an  elastomer1^,  by  using  an  elastomer  of  interpenetrating 
networks^, 55f  by  swelling  an  elastomer  with  a  solvent^6!,  and  by  spraying  charge  on  an 
electrode-free  elastomer  [57]. 

These  experimental  advances  have  prompted  a  theoretical  question:  What  is  the 
fundamental  limit  of  deformation  that  can  be  induced  by  voltage?  One  can  easily  increase  the 
length  of  a  rubber  band  several  times  by  using  a  mechanical  force.  Why  is  it  difficult  to  do  so  by 
using  a  voltage?  The  difficulty  has  to  do  with  two  modes  of  failure  associated  with  apply  a 
voltage:  electrical  breakdown  and  electromechanical  instability.  For  a  stiff  dielectric  such  as  a 
ceramic  or  a  glassy  polymer,  voltage-induced  deformation  is  limited  by  electrical  breakdown, 
when  the  voltage  mobilizes  charged  species  in  the  dielectric  to  produce  a  path  of  electrical 
conduction. 

For  a  compliant  dielectric  such  as  an  elastomer,  the  voltage-induced  deformation  is 
often  limited  by  electromechanical  instability.  It  was  Stark  and  Gartonfss]  who  first  reported 
that  the  breakdown  fields  of  polymers  reduced  when  the  polymers  became  soft  at  elevated 
temperatures.  The  phenomenon  is  understood  as  follows.  The  electric  voltage  is  applied 
between  the  electrodes  on  the  top  and  the  bottom  surfaces  of  a  thin  layer  of  a  polymer.  As  the 
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electric  field  increases,  the  polymer  thins  down,  so  that  the  same  voltage  will  induce  an  even 
higher  electric  field.  This  positive  feedback  results  in  a  mode  of  instability,  known  as 
electromechanical  instability  or  pull-in  instability,  which  causes  the  polymer  to  reduce  the 
thickness  drastically,  often  leading  to  electrical  breakdown.  Electromechanical  instability  has 
been  recognized  as  a  failure  mode  of  the  insulators  for  power  transmission  cables. 

Electromechanical  instability  is  sensitive  to  the  stress-stretch  behavior  of  the 
elastomer^].  Fig.  12a  sketches  a  dielectric  membrane  pulled  by  biaxial  stresses  a .  The  length 
of  the  membrane  in  any  direction  in  the  plane  is  stretched  by  a  ratio  A .  As  will  become  clear,  to 
attain  a  large  voltage-induced  stretch,  the  dielectric  should  have  a  stress-stretch  curve  <j(a)  of 
the  following  desirable  features^91:  (a)  The  dielectric  is  compliant  at  small  stretches,  and  (b)  the 
dielectric  stiffens  steeply  at  modest  stretches.  That  is,  the  limiting  stretch,  Aljm ,  should  not  be 
excessive.  Also  sketched  are  several  designs  of  materials  that  exhibit  the  stress-stretch  curve  of 
the  desirable  form.  Many  biological  tissues,  such  as  skins  and  vascular  walls,  deform  readily, 
but  avert  excessive  deformation.  Fig.  12b  sketches  a  design  of  such  a  tissue,  consisting  of  stiffer 
fibers  in  a  compliant  matrix.  At  small  stretches,  the  fibers  are  loose,  and  the  tissue  is  compliant. 
At  large  stretches,  the  fibers  are  taut,  and  the  tissue  stiffens  steeply.  As  another  example,  Fig. 
12c  sketches  a  network  of  polymers  with  folded  domains.  The  domains  unfold  when  the 
network  is  pulled,  giving  rise  to  substantial  deformation.  After  all  the  domains  unfold,  the 
network  stiffens  steeply. 

Consider  a  synthetic  elastomer,  i.e.,  a  network  of  polymer  chains.  When  the  individual 
chains  are  short,  the  initial  modulus  of  the  elastomer  is  large  and  the  limiting  stretch  AVim  is 
small.  When  the  individual  chains  are  long,  the  initial  modulus  of  the  elastomer  is  small  and  the 
limiting  stretch  A]im  is  large.  Consequently,  it  is  difficult  to  achieve  the  stress-stretch  curve  of 
the  desirable  form  by  adjusting  the  density  of  crosslinks  alone.  The  stress-stretch  curve, 
however,  can  be  shaped  into  the  desirable  form  in  several  ways.  For  example,  the  widely  used 
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dielectric  elastomer,  VHB,  is  a  network  of  polymers  with  side  chains  (Fig.  I2d).  The  side  chains 
fill  the  space  around  the  networked  chains.  The  motion  of  the  networked  chains  is  lubricated, 
lowering  the  glass  transition  temperature.  Also  the  density  of  the  networked  chains  is  reduced, 
lowering  the  stiffness  of  the  elastomer  when  the  stretch  is  small.  While  the  side  chains  do  not 
change  the  contour  length  of  the  networked  chains,  the  side  chains  pull  the  networked  chains 
towards  their  full  contour  length  even  when  the  elastomer  is  not  loaded.  Once  loaded,  the 
elastomer  may  stiffen  sharply,  averting  electromechanical  instability.  Similar  behavior  is 
expected  for  a  network  swollen  with  a  solvent  (Fig.  I2e).  The  stress-stretch  curve  can  also  be 
shaped  into  the  desirable  form  by  prestretch1!1,  or  by  using  interpenetrating  networks1 [54,55]. 

As  illustrated  in  Fig.  13a,  when  a  membrane  of  an  elastomer,  thickness  H  in  the 
undeformed  state,  is  subject  to  a  voltage  ® ,  the  membrane  is  stretched  by  A  in  both  directions 
in  the  plane,  the  thickness  of  the  membrane  reduces  to  HA  2 ,  and  the  electric  field  in  the 
membrane  is  E  -  T2®  /  H .  The  membrane  is  taken  to  be  incompressible.  The  actuation  can  be 
described  by  the  Maxwell  stress.  A  combination  of  the  above  considerations  relates  the  voltage 
to  the  stretch: 

®  =  /M2VexU)/i- .  (70) 

This  voltage-stretch  relation  is  sketched  in  Fig.  13a.  Even  though  the  stress-stretch  curve  <j(a) 
is  monotonic,  the  voltage-stretch  curve  ®(/l)  is  usually  not[39l.  At  a  small  stretch  (A  ~  1),  the 
rising  <j(a)  dominates,  and  the  voltage  increases  with  the  stretch.  At  an  intermediate  stretch, 

the  factor  X*  due  to  thinning  of  the  membrane  becomes  important,  and  the  voltage  falls  as  the 
stretch  increases.  As  the  elastomer  approaches  the  limiting  stretch  Alim ,  the  steep  rise  of  ct(a) 
prevails,  and  the  voltage  rises  again.  The  shape  of  the  voltage-stretch  curve  ®(t)  indicates  a 
snap-through  electromechanical  instability^391.  The  instability  can  cause  some  regions  of  the 
elastomer  to  thin  down  more  than  others. 
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The  local  maximum  voltage  represents  a  critical  condition,  which  can  be  estimated  as 
follows.  Under  the  equal -biaxial  stresses,  Hooke’s  law  takes  the  form  where  // 

is  the  shear  modulus.  Inserting  this  expression  into  (70),  and  maximizing  the  function  o(k) , 

we  find  local  maximum  voltage  ®c  «  0.8 oil -yjju/s  and  the  critical  4  =4/3  =  1.3  3.  The  critical 
values  vary  somewhat  with  the  stress-stretch  relation.  For  example,  for  the  neo-Hookean  model, 
<t(t)  =  /j(a2  — /T4) ,  the  maximum  voltage  is  ®c « 0.6  gHyfjuJs  and  the  critical  stretch  is 
Xc  =  21/3  « 1.2  6.  This  electromechanical  instability  has  been  analyzed  systematically  by  using  the 
Hessian  [6<>64]. 

Before  a  voltage  is  applied,  an  elastomer  may  be  prestretched  to  A  by  a  mechanical  force, 
and  then  fixed  by  rigid  electrodes.  Subsequently,  when  the  voltage  is  applied,  the  elastomer  will 
not  deform  further.  The  measured  voltage  at  failure  is  taken  to  be  the  electrical  breakdown 
voltage.  Experiments  indicate  that  the  breakdown  voltage  is  a  monotonically  decreasing 
function  of  the  prestretch,  ®,,  (7)  [30,41]. 

According  to  where  the  curves  ®U)  and  ® ,,(/!.  )  intersect,  we  distinguish  three  types  of 
dielectric  elastomers [59].  A  type  I  dielectric  suffers  electrical  breakdown  prior  to 
electromechanical  instability,  and  is  capable  of  small  deformation  of  actuation,  Fig.  13b.  A  type 
II  dielectric  reaches  the  peak  of  the  ®(/l)  curve,  and  thins  down  excessively,  leading  to  electrical 
breakdown,  Fig.  13c.  The  dielectric  is  recorded  to  fail  at  the  peak  of  ®(/l) ,  which  can  be  much 
below  the  breakdown  voltage  ®  B .  The  deformation  of  actuation  is  limited  by  the  stretch  at 
which  the  voltage  reaches  the  peak.  A  type  III  dielectric  eliminates  or  survives 
electromechanical  instability,  reaches  a  stable  state  before  the  electrical  breakdown,  and  attains 
a  large  deformation  of  actuation,  Fig.  13d. 

A  new  experimental  manifestation  of  the  electromechanical  instability  has  been 
reported  recently[3°].  Under  certain  conditions,  an  electric  voltage  can  deform  a  layer  of  a 
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dielectric  elastomer  into  a  mixture  of  two  regions,  one  being  flat  and  the  other  wrinkled  (Fig.  14). 
In  the  experiment,  the  electrodes  on  the  top  and  the  bottom  surfaces  of  the  dielectric  layer  were 
made  of  conducting  grease,  which  applied  a  uniform  electric  potential  to  the  elastomer  without 
constraining  its  deformation.  This  observation  has  been  interpreted  as  the  coexistence  of  two 
states[39]. 


VI.  INHOMOGENEOUS  FIELDS 

Studies  of  inhomogeneous  fields  of  coupled  large  deformation  and  electric  potential  date 
back  to  classic  works  of  Toupin[66],  Eringen[67]  and  Tiersten[68f  These  works  have  been 
reexamined  recently  for  applications  to  dielectric  elastomers  [9'13].  Here  we  summarize  basic 
ideas,  following  the  presentation  of  Ref.  [12]. 


6.1.  Condition  of  Thermodynamic  Equilibrium 

A  body  of  an  elastic  dielectric  is  represented  by  a  field  of  material  particles.  Each 
material  particle  is  named  after  the  coordinate  X  of  its  place  when  the  body  is  in  a  reference 
state.  In  the  current  state,  at  time  t,  the  particle  X  moves  to  a  place  with  coordinate  x.  The 
function 

x  =  x(X,f)  (71) 


describes  the  history  of  the  deformation  of  the  body.  Define  the  deformation  gradient  F  as 


ar  ,.(x,0 

riK  — 


dX 


(72) 


K 


The  deformation  gradient  generalizes  the  notion  of  the  stretches. 

In  the  current  state  at  time  t,  the  electric  potential  at  particle  X  is  denoted  as 

cp  =  o(x,t). 


(73) 


The  gradient  of  the  electric  potential  defines  the  nominal  electric  field  E  ,  namely, 
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(74) 


E 


K 


50  (X,f) 

sxK 


The  negative  sign  in  (74)  follows  the  convention  that  the  electric  field  vector  points  in  the 
direction  from  a  particle  of  a  high  voltage  to  a  particle  of  a  low  voltage. 

Motivated  by  (22),  we  write  the  variation  of  the  nominal  Helmholtz  free  energy,  SW  ,  in 
the  form: 


SW  =  sjKSFiK 


+ eksdk  , 


(75) 


where  SFiK  is  a  small  change  in  the  deformation  gradient,  and  SDK  is  a  small  change  in  the 
nominal  electric  field.  Equation  (75)  defines  the  nominal  stress  s  as  a  tensor  work-conjugate  to 
the  deformation  gradient  F,  and  the  nominal  electric  displacement  D  as  a  vector  work- 
conjugate  to  the  nominal  electric  field  E  . 

Inspecting  (72)  and  (74),  we  wish  to  use  the  deformation  gradient  and  the  nominal 
electric  field  as  the  independent  variables.  Introducing  a  new  quantity  W  by 

W  =  W-EkDk.  (76) 


The  quantity  W  may  be  called  the  electrical  Gibbs  free  energy.  A  combination  of  (75)  and  (76) 
gives 

SW  =  s!KdEK-DKdEK.  (77) 

We  may  call  the  quantity  DK8EK  the  complementary  electrical  work. 

A  material  model  is  prescribed  by  a  function  W  -  w(f,e).  When  the  body  undergoes  a 
rigid  body  motion,  the  free  energy  is  invariant.  Consequently,  the  function  depends  on  the 
deformation  gradient  through  the  Green  deformation  tensor,  CK  =  FiKFiL .  Associated  with 


small  changes  SFiK  and  5EA. ,  the  electrical  Gibbs  free  energy  changes  by 


-  dwi F  e)  dwi F  E  ~ 

SW  =  \  ’  >  8F.k+  '  5Ek. 

dFiv  5E, 


(78) 
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On  each  material  element  of  volume  dV(x) ,  we  prescribe  mass  p{X)dV  ,  force 
B(X,f)dV  and  charge  q(X,t)dV .  The  effect  of  inertia  may  be  represented  by  adding  to  the 
force  the  inertial  force,  so  that  the  net  force  on  the  element  of  volume  is  (  /!,  -  pdf  x- /  dt'z  ]de,d  V . 
On  each  material  element  of  interface  dA(x) ,  we  prescribe  force  T(X,Z  )dA  and  charge 
co(X.,t)dA. 

Let  3c(  =A,.(x)  be  a  field  of  virtual  displacement  of  the  body— that  is,  every  material 
particle  X  moves  independently.  The  virtual  displacement  A,.(X)  is  unrelated  to  actual 
displacement  x;(X,t).  Associated  with  the  field  of  virtual  displacement,  the  forces  do  virtual 

work  |(s,  - pd*xj dt^yk^V +^Ti&cidA.  Similarly,  let  dP  =  //(X)  be  a  field  of  virtual  electric 
potential  of  the  body.  Associated  with  the  field  of  electric  potential,  the  charges  do  virtual 
complementary  work  I  gdPdF  + 1  coS&dA  .  The  virtual  deformation  gradient  is 


3FiK=dAi{x)/ dXK  ,  the  virtual  nominal  electric  field  is  5EK=-8q{X)/ dXK  ,  and  the  virtual 

r  /v  a. 

change  in  the  electrical  Gibbs  free  energy  is  SWdV ,  where  SW  is  given  by  (78).  When  the 

body  is  in  thermodynamic  equilibrium,  the  change  in  the  electrical  Gibbs  free  energy  equals  the 
mechanical  work  minus  the  complementary  electrical  work: 


p\  2  \ 

JdHW  =  J  &;dF  +  J T)<3vr;dA- j" qcXPdF- j" cod&dA . 


dt 


(79) 


This  condition  of  thermodynamic  equilibrium  has  the  similar  physical  content  as  (1)  and  (21), 
and  holds  for  arbitrary  and  independent  variations  dx  and  dP . 

Once  the  loads  and  the  electrical  Gibbs  free-energy  function  w[f,  e)  are  prescribed,  the 
variational  statement  (79),  along  with  the  definitions  (72)  and  (74),  is  the  basis  for  the  finite 
element  method,  determining  the  field  of  deformation  x(X,f)  and  the  field  of  electric  potential 
<p(X,f)  simultaneously.  Several  implementations  of  the  finite  element  method  have  been 
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reported16^69'7'!,  but  few  practical  examples  are  available,  especially  transducers  approach 
electromechanical  instability.  Significant  effort  is  needed  to  develop  the  finite  element  method, 
and  to  apply  the  method  to  analyze  phenomena  and  devices. 


6.2.  Differential  Equations 

The  variational  statement  of  thermodynamic  equilibrium  also  leads  to  partial  differential 
equations.  These  equations  are  listed  in  this  subsection. 

A  comparison  of  (77)  and  (78)  gives  that 


s 


iK 


dw( F,E 

dF 

U1  iK 


D 


K 


8w{f,E 

8Ek 


(80) 


(81) 


Once  the  electrical  free-energy  function  Iv(f,e)  is  prescribed,  (80)  and  (81)  constitute  the 
equations  of  states. 

Inserting  (72),  (74)  and  (77)  into  the  condition  of  thermodynamic  equilibrium  (79),  and 
recalling  that  the  condition  holds  for  arbitrary  and  independent  variations  in  dx  and  & f> ,  we 
obtain  that 


dsJXj) 

dXK 


+  B,(X,t)=p(X) 


d%{X,t) 


dt2 


(82) 


in  the  volume 


on  an  interface, 


=Tt 


dPK{X,t ) 

dXK 


q(x,f) 


(83) 


(84) 


in  the  volume,  and 


(D+K-DK)NK=co{X,t) 


(85) 
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on  the  interfaces.  Equations  (82)  and  (83)  reproduce  the  equations  for  momentum  balance, 
and  (84)  and  (85)  reproduce  Gauss’s  law  of  electrostatics. 

Equations  (71M74)  and  (80M85)  are  governing  equations  to  determine  the  field  of 
deformation  x(X,t)  and  the  field  of  electric  potential  ®(x,f)  simultaneously,  once  the  loads 

and  the  free-energy  function  1e(f,e)  are  prescribed.  These  partial  differential  equations  have 
been  used  to  solve  boundary-value  problems  involving  coupled  large  deformation  electric 
potential1^2-?6].  observe  that  the  equations  of  mechanics,  (71),  (72),  (82)  and  (83),  decouple 
from  those  of  electrostatics,  (73),  (74),  (84)  and  (85).  The  only  coupling  between  mechanics 
and  electrostatics  arises  from  the  material  model,  (80)  and  (81). 


6.3.  True  Quantities 

The  true  stress  cr.  relates  to  the  nominal  stress  by 

Fjk 

detF  S'7'  ‘ 

The  true  electric  displacement  Dj  relates  to  the  nominal  electric  displacement  as 

F  ~ 

D;=—^Dk. 

'  detF  * 

The  true  electric  field  E,  relates  to  the  nominal  electric  field  as 


(86) 


(87) 


E^HiKEK,  (88) 

where  HiK  is  the  inverse  of  the  deformation  gradient,  namely,  HiKFjL  = s kl  and  HiKFjK  =  S... 

The  true  quantities  are  functions  of  x  and  t,  and  satisfy  the  familiar  partial  differential  equations 
in  mechanics  and  electrostatics. 


6.4.  Ideal  Dielectric  Elastomers 

Toupin[66]  noted  that,  for  an  isotropic  elastic  dielectric,  the  free-energy  density  is  a 
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function  of  six  invariants  of  the  deformation  gradient  tensor  and  the  electric  field  vector. 
Function  of  this  complexity  is  unavailable  for  any  real  material.  Several  further  considerations 
may  reduce  the  complexity  of  the  free-energy  function  somewhat,  but  are  still  far  from  being 
useful  in  practice[12]. 

We  next  describe  ideal  dielectric  elastomers— a  material  model  nearly  exclusively  used  in 
the  literature.  As  discussed  before  in  connection  with  Fig.  9,  for  an  ideal  dielectric  elastomer, 
the  dielectric  behavior  is  the  same  as  that  of  a  liquid— that  is,  the  dielectric  behavior  is 
unaffected  by  deformation1^].  As  a  simplest  model  of  a  dielectric  liquid,  assume  that  the  true 
electric  displacement  Dm  is  linear  in  the  true  electric  field  Em  : 

Dm=sEm  (89) 

The  permittivity  s  is  taken  to  be  independent  of  deformation. 

Using  (87)  and  (88),  we  express  (89)  in  terms  of  the  nominal  fields: 

DN  =  £E,  HmNHrnI'  detF .  (90) 

Inserting  (90)  into  (81)  and  integrating,  we  obtain  the  nominal  density  of  the  electrical  Gibbs 
free  energy: 

M F,e)=  detF .  (91) 

Here  WS(F)  is  the  free  energy  associated  with  the  strength  of  the  elastomer,  which  may  be 
selected  from  a  large  menu  in  the  theory  of  elasticity.  While  the  elastomer  is  nearly 
incompressible,  in  the  finite  element  method,  it  is  convenient  to  allow  the  material  to  be 
compressible  with  a  large  bulk  modulus. 

Insert  (91)  into  (80),  and  recall  mathematical  identities  oHmN  / 8 FiK  =  -HmKHjN  and 
o  detF  /  8FiK  =  HiK  detF .  We  obtain  that 

5W(f)  ~  ~  (  1  A 

siK  =  +  sEnEl  HiNHmKHnJ  -  -  HmNHmLHiK  detF .  (92) 

SFiK  K  2 
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This  equation  of  state  expresses  the  nominal  stress  as  a  function  of  the  deformation  gradient 
and  the  nominal  electric  field. 

A  combination  of  (86)  and  (92)  gives 


CT  fk 

v  detF  dF„ 


EE  ■ - E,„E,„d:: 


(93) 


This  equation  expresses  the  true  stress  as  a  function  of  the  deformation  gradient  and  the  true 
electric  field.  The  contribution  due  to  the  deformation  gradient  results  of  the  stretching  of  the 
elastomer,  and  is  the  same  as  that  in  the  theory  of  elasticity.  The  contribution  due  to  the  electric 
field  is  identical  to  that  derived  by  Maxwell[40].  When  cr  =  o ,  (93)  balances  elasticity  and 


electrostatics,  and  determines  the  voltage-induced  deformation.  As  commented  before,  the 
Maxwell  stress  correctly  accounts  for  the  voltage-induced  deformation  only  when  the  dielectric 
behavior  is  liquid-like,  an  idealization  works  well  with  elastomers,  but  not  for  any  other  solid 
dielectrics. 

Equations  (89)  and  (93)  constitute  the  equations  of  state,  for  an  ideal  dielectric 
elastomer,  in  terms  of  the  true  quantities.  The  equations  of  state  exhibit  one-way  coupling:  the 
deformation  does  not  affect  the  dielectric  behavior,  but  the  electric  field  contributes  to  the 
stress-stretch  relation.  As  noted  before,  the  partial  differential  equations  of  mechanics 
decouples  from  those  of  electrostatics.  One  may  solve  electrostatic  boundary-value  problems  in 
terms  of  the  true  fields,  and  then  add  the  Maxwell  stress  in  solving  the  elastic  field.  Of  course, 
the  deformation  will  change  the  shape  of  the  boundaries  of  the  body.  This  change  must  be 
included  in  solving  the  electrostatic  problems.  The  one-way  coupling  may  not  bring  any 
advantage  after  all. 

The  model  of  ideal  dielectric  elastomers  can  be  generalized  to  account  for  nonlinear 
dielectric  behavior,  by  replacing  (89)  with  a  nonlinear  relation  between  the  electric  field  and 
electric  displacement1'21.  Nonlinear  dielectric  behavior  may  be  significant  at  high  electric  fields. 
Furthermore,  the  model  of  ideal  dielectric  elastomers  can  be  modified  to  include  dissipative 
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processes,  such  as  viscoelasticity,  dielectric  relaxation,  and  electrical  conductions1^0]. 


VII.  CONCLUDING  REMARKS 

A  large  number  of  examples  in  biology  demonstrate  that  deformation  of  soft  materials 
connect  diverse  stimuli  to  many  functions  essential  to  the  life.  An  exciting  field  of  engineering  is 
emerging  that  uses  soft  active  materials  to  create  soft  machines.  To  participate  in  advancing  the 
field  of  soft  active  materials  and  soft  machines  effectively,  mechanicians  must  retool  our 
laboratories  and  our  software,  as  well  as  adapt  our  theories.  While  theories  are  being  developed 
for  diverse  soft  active  materials,  this  review  focuses  on  one  class  of  soft  active  materials: 
dielectric  elastomers.  This  focus  allows  us  to  review  the  theory  in  some  depth,  within  the 
framework  of  nonlinear  continuum  mechanics  and  nonequilibrium  thermodynamics,  while 
motivating  the  theory  by  empirical  observations,  molecular  pictures  and  applications.  It  is 
hoped  that  the  theory  will  be  used  to  develop  software,  study  intriguing  phenomena  of 
electromechanical  coupling,  and  aid  the  creation  of  soft  active  materials  and  soft  machines. 
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Stimuli  cause  deformation 


Deformation  provide  functions 


Fig.i.  The  environment  affects  a  material  through  diverse  stimuli,  such  as  a  force,  an  electric 
field,  a  charge  in  pH,  and  a  change  in  temperature.  In  response  to  a  stimulus,  a  soft  active 
material  (SAM)  deforms.  The  deformation  provides  a  function,  such  as  change  in  color  and 
change  in  flow  rate. 
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Reference  state 


Fig.  2.  A  dielectric  elastomer  in  the  reference  state  and  in  a  current  state. 
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electrode 


Fig.3.  A  transducer  consists  of  a  dielectric  separating  two  electrodes.  The  transducer  is  subject 
to  a  force,  represented  by  a  weight  P.  The  two  electrodes  are  connected  through  a  conducting 
wire  to  a  battery  of  voltage  ®  .  The  weight  moves  by  distance  l,  and  the  battery  pumps  charge  Q 
from  one  electrode  to  the  other. 
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o  =  o1 


Fig.  4.  In  a  plane  with  force  and  displacement  as  coordinates,  a  point  represents  a  state  of  a 
transducer.  A  curve  of  a  constant  voltage  is  the  force-displacement  curve  measured  when  the 
transducer  is  subject  to  a  constant  voltage. 
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Fig.  5.  In  a  plane  with  voltage  and  charge  as  coordinates,  a  point  represents  a  state  of  a 
transducer.  A  use  of  the  transducer  typically  involves  a  cyclic  change  of  the  state.  The  rectangle 
represents  a  cycle  involving  two  levels  of  voltage  and  two  values  of  charge. 
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FIG.  6.  A  state  of  a  dielectric  membrane  is  represented  by  a  point  in  the  charge-voltage  plane. 
The  coordinates  are  given  in  dimensionless  forms.  Plotted  are  curves  representing  various 
modes  of  failure:  electrical  breakdown  (EB),  electromechanical  instability  (EMI),  loss  of  tension 
(s  =  o),  and  rupture  by  stretch  ( X  =  XR ).  These  curves  bound  the  region  of  allowable  states  of 
the  transducer.  A  cycle  involving  two  levels  of  voltage  and  two  values  of  charge  is  represented 
by  dotted  lines. 
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Fig.  7.  A  parallel-plate  capacitor  consists  of  two  electrodes  separated  by  a  thin  gap  of  a  vacuum. 
When  a  voltage  is  applied,  the  two  electrodes  attract  each  other.  The  electrostatic  attraction  is 
balanced  by  applying  a  force. 
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Fig.  8.  A  parallel-plate  capacitor  is  partly  in  the  air  and  partly  in  a  liquid.  When  a  voltage  is 
applied,  the  liquid  rises.  As  indicated  by  the  free-body  diagram  on  the  right,  the  rise  of  the 
liquid  is  due  to  the  balance  between  the  Maxwell  stress  and  the  weight  of  the  liquid. 
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Fig.  9.  An  elastomer  is  a  three  dimensional  network  of  long  and  flexible  polymer  chains.  Each 
polymer  chain  consists  of  a  large  number  of  monomers. 
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triaxial  stresses  uniaxial  biaxial  stresses 


Fig.io.  A  dielectric  in  three  states  of  stresses. 
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(a) 


(c) 


+++++++++ 


(b) 


Fig.  11.  Consider  a  dielectric  that  is  nonpolar  in  the  absence  of  applied  voltage  (a).  Subject  to 
voltage,  some  dielectrics  become  thinner  (b),  but  other  dielectrics  become  thicker  (c). 
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FIG.  12.  (a)  Stress-stretch  curve  of  a  membrane  under  biaxial  stresses,  (b)  Fibers  embedded  in 
a  compliant  matrix,  (c)  A  network  of  polymers  with  folded  domains,  (d)  A  network  of  polymers 
with  side  chains,  (e)  A  network  of  polymers  swollen  with  a  solvent. 
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FIG.  13.  (a)  A  membrane  of  a  dielectric  elastomer  subject  to  a  voltage  reduces  thickness  and 
expands  area.  The  voltage-stretch  curve  is  typically  not  monotonic,  (b-d)  Three  types  behavior 
are  distinguished,  depending  on  where  the  two  curves  0(1)  and  O ,,  (/.)  intersect. 
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Fig.  14.  An  experimental  observation  of  electromechanical  instability  (courtesy  of  JS  Plante  and 
S  Dubowsky).  A  layer  of  a  dielectric  elastomer,  coated  with  conductive  grease  on  top  and  bottom 
faces,  is  pre-stretched  using  a  frame.  An  electric  voltage  is  applied  between  the  two  electrodes. 
The  layer  deforms  into  a  mixture  of  two  regions,  one  being  flat  and  the  other  wrinkled. 
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